Heat shock protein 90 facilitates SARS-CoV-2 structural protein-mediated virion assembly and promotes virus-induced pyroptosis.

Inhibition of heat shock protein 90 (Hsp90), a prominent molecular chaperone, effectively limits severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection but little is known about any interaction between Hsp90 and SARS-CoV-2 proteins. Here, we systematically analyzed the effects of the chaperone isoforms Hsp90α and Hsp90ß on individual SARS-CoV-2 viral proteins. Five SARS-CoV-2 proteins, namely nucleocapsid (N), membrane (M), and accessory proteins Orf3, Orf7a, and Orf7b were found to be novel clients of Hsp90ß in particular. Pharmacological inhibition of Hsp90 with 17-DMAG results in N protein proteasome-dependent degradation. Hsp90 depletion-induced N protein degradation is independent of CHIP, a ubiquitin E3 ligase previously identified for Hsp90 client proteins, but alleviated by FBXO10, an E3 ligase identified by subsequent siRNA screening. We also provide evidence that Hsp90 depletion may suppress SARS-CoV-2 assembly partially through induced M or N degradation. Additionally, we found that GSDMD-mediated pyroptotic cell death triggered by SARS-CoV-2 was mitigated by inhibition of Hsp90. These findings collectively highlight a beneficial role for targeting of Hsp90 during SARS-CoV-2 infection, directly inhibiting virion production and reducing inflammatory injury by preventing the pyroptosis that contributes to severe SARS-CoV-2 disease.

Distinct Conformations of SARS-CoV-2 Omicron Spike Protein and Its Interaction with ACE2 and Antibody.

Since November 2021, Omicron has been the dominant severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant that causes the coronavirus disease 2019 (COVID-19) and has continuously impacted human health. Omicron sublineages are still increasing and cause increased transmission and infection rates. The additional 15 mutations on the receptor binding domain (RBD) of Omicron spike proteins change the protein conformation, enabling the Omicron variant to evade neutralizing antibodies. For this reason, many efforts have been made to design new antigenic variants to induce effective antibodies in SARS-CoV-2 vaccine development. However, understanding the different states of Omicron spike proteins with and without external molecules has not yet been addressed. In this review, we analyze the structures of the spike protein in the presence and absence of angiotensin-converting enzyme 2 (ACE2) and antibodies. Compared to previously determined structures for the wildtype spike protein and other variants such as alpha, beta, delta, and gamma, the Omicron spike protein adopts a partially open form. The open-form spike protein with one RBD up is dominant, followed by the open-form spike protein with two RBD up, and the closed-form spike protein with the RBD down. It is suggested that the competition between antibodies and ACE2 induces interactions between adjacent RBDs of the spike protein, which lead to a partially open form of the Omicron spike protein. The comprehensive structural information of Omicron spike proteins could be helpful for the efficient design of vaccines against the Omicron variant.

Imprinted SARS-CoV-2 humoral immunity induces convergent Omicron RBD evolution.

Continuous evolution of Omicron has led to a rapid and simultaneous emergence of numerous variants that display growth advantages over BA.5 (ref. 1). Despite their divergent evolutionary courses, mutations on their receptor-binding domain (RBD) converge on several hotspots. The driving force and destination of such sudden convergent evolution and its effect on humoral immunity remain unclear. Here we demonstrate that these convergent mutations can cause evasion of neutralizing antibody drugs and convalescent plasma, including those from BA.5 breakthrough infection, while maintaining sufficient ACE2-binding capability. BQ.1.1.10 (BQ.1.1 + Y144del), BA.4.6.3, XBB and CH.1.1 are the most antibody-evasive strains tested. To delineate the origin of the convergent evolution, we determined the escape mutation profiles and neutralization activity of monoclonal antibodies isolated from individuals who had BA.2 and BA.5 breakthrough infections2,3. Owing to humoral immune imprinting, BA.2 and especially BA.5 breakthrough infection reduced the diversity of the neutralizing antibody binding sites and increased proportions of non-neutralizing antibody clones, which, in turn, focused humoral immune pressure and promoted convergent evolution in the RBD. Moreover, we show that the convergent RBD mutations could be accurately inferred by deep mutational scanning profiles4,5, and the evolution trends of BA.2.75 and BA.5 subvariants could be well foreseen through constructed convergent pseudovirus mutants. These results suggest that current herd immunity and BA.5 vaccine boosters may not efficiently prevent the infection of Omicron convergent variants.

Thermal Titration Molecular Dynamics (TTMD): Not Your Usual Post-Docking Refinement.

Molecular docking is one of the most widely used computational approaches in the field of rational drug design, thanks to its favorable balance between the rapidity of execution and the accuracy of provided results. Although very efficient in exploring the conformational degrees of freedom available to the ligand, docking programs can sometimes suffer from inaccurate scoring and ranking of generated poses. To address this issue, several post-docking filters and refinement protocols have been proposed throughout the years, including pharmacophore models and molecular dynamics simulations. In this work, we present the first application of Thermal Titration Molecular Dynamics (TTMD), a recently developed method for the qualitative estimation of protein-ligand unbinding kinetics, to the refinement of docking results. TTMD evaluates the conservation of the native binding mode throughout a series of molecular dynamics simulations performed at progressively increasing temperatures with a scoring function based on protein-ligand interaction fingerprints. The protocol was successfully applied to retrieve the native-like binding pose among a set of decoy poses of drug-like ligands generated on four different pharmaceutically relevant biological targets, including casein kinase 1δ, casein kinase 2, pyruvate dehydrogenase kinase 2, and SARS-CoV-2 main protease.

The Cytoplasmic Domain of the SARS-CoV-2 Envelope Protein Assembles into a ß-Sheet Bundle in Lipid Bilayers.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) envelope (E) protein forms a pentameric ion channel in the lipid membrane of the endoplasmic reticulum Golgi intermediate compartment (ERGIC) of the infected cell. The cytoplasmic domain of E interacts with host proteins to cause virus pathogenicity and may also mediate virus assembly and budding. To understand the structural basis of these functions, here we investigate the conformation and dynamics of an E protein construct (residues 8-65) that encompasses the transmembrane domain and the majority of the cytoplasmic domain using solid-state NMR. 13C and 15N chemical shifts indicate that the cytoplasmic domain adopts a ß-sheet-rich conformation that contains three ß-strands separated by turns. The five subunits associate into an umbrella-shaped bundle that is attached to the transmembrane helices by a disordered loop. Water-edited NMR spectra indicate that the third ß-strand at the C terminus of the protein is well hydrated, indicating that it is at the surface of the ß-bundle. The structure of the cytoplasmic domain cannot be uniquely determined from the inter-residue correlations obtained here due to ambiguities in distinguishing intermolecular and intramolecular contacts for a compact pentameric assembly of this small domain. Instead, we present four structural topologies that are consistent with the measured inter-residue contacts. These data indicate that the cytoplasmic domain of the SARS-CoV-2 E protein has a strong propensity to adopt ß-sheet conformations when the protein is present at high concentrations in lipid bilayers. The equilibrium between the ß-strand conformation and the previously reported α-helical conformation may underlie the multiple functions of E in the host cell and in the virion.

Interactions of SARS-CoV-2 and MERS-CoV fusion peptides measured using single-molecule force methods.

We address the challenge of understanding how hydrophobic interactions are encoded by fusion peptide (FP) sequences within coronavirus (CoV) spike proteins. Within the FPs of severe acute respiratory syndrome CoV 2 and Middle East respiratory syndrome CoV (MERS-CoV), a largely conserved peptide sequence called FP1 (SFIEDLLFNK and SAIEDLLFDK in SARS-2 and MERS, respectively) has been proposed to play a key role in encoding hydrophobic interactions that drive viral-host cell membrane fusion. Although a non-polar triad (Leu-Leu-Phe (LLF)) is common to both FP1 sequences, and thought to dominate the encoding of hydrophobic interactions, FP1 from SARS-2 and MERS differ in two residues (Phe 2 versus Ala 2 and Asn 9 versus Asp 9, respectively). Here we explore whether single-molecule force measurements can quantify hydrophobic interactions encoded by FP1 sequences, and then ask whether sequence variations between FP1 from SARS-2 and MERS lead to significant differences in hydrophobic interactions. We find that both SARS-2 and MERS wild-type FP1 generate measurable hydrophobic interactions at the single-molecule level, but that SARS-2 FP1 encodes a substantially stronger hydrophobic interaction than its MERS counterpart (1.91 ± 0.03 nN versus 0.68 ± 0.03 nN, respectively). By performing force measurements with FP1 sequences with single amino acid substitutions, we determine that a single-residue mutation (Phe 2 versus Ala 2) causes the almost threefold difference in the hydrophobic interaction strength generated by the FP1 of SARS-2 versus MERS, despite the presence of LLF in both sequences. Infrared spectroscopy and circular dichroism measurements support the proposal that the outsized influence of Phe 2 versus Ala 2 on the hydrophobic interaction arises from variation in the secondary structure adopted by FP1. Overall, these insights reveal how single-residue diversity in viral FPs, including FP1 of SARS-CoV-2 and MERS-CoV, can lead to substantial changes in intermolecular interactions proposed to play a key role in viral fusion, and hint at strategies for regulating hydrophobic interactions of peptides in a range of contexts.

Variations within the Glycan Shield of SARS-CoV-2 Impact Viral Spike Dynamics.

The emergence of SARS-CoV-2 variants alters the efficacy of existing immunity, whether arisen naturally or through vaccination. Understanding the structure of the viral spike assists in determining the impact of mutations on the antigenic surface. One class of mutation impacts glycosylation attachment sites, which have the capacity to influence the antigenic structure beyond the immediate site of attachment. Here, we compare the site-specific glycosylation of recombinant viral spike mimetics of B.1.351 (Beta), P.1 (Gamma), B.1.617.2 (Delta), B.1.1.529 (Omicron). The P.1 strain exhibits two additional N-linked glycan sites compared to the other variants analyzed and we investigate the impact of these glycans by molecular dynamics. The acquired N188 site is shown to exhibit very limited glycan maturation, consistent with limited enzyme accessibility. Structural modeling and molecular dynamics reveal that N188 is located within a cavity by the receptor binding domain, which influences the dynamics of these attachment domains. These observations suggest a mechanism whereby mutations affecting viral glycosylation sites have a structural impact across the protein surface.

Potential targets of severe acute respiratory syndrome coronavirus 2 of clinical drug fluvoxamine: Docking and molecular dynamics studies to elucidate viral action.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has continued evolving for survival and adaptation by mutating itself into different variants of concern, including omicron. Several studies and clinical trials found fluvoxamine, an Food and Drug Administration-approved antidepressant drug, to be effective at preventing mild coronavirus disease 2019 (COVID-19) from progressing to severe diseases. However, the mechanism of fluvoxamine's direct antiviral action against COVID-19 is still unknown. Fluvoxamine was docked with 11 SARS-CoV-2 targets and subjected to stability, conformational changes, and binding free energy analyses to explore its mode of action. Of the targets, nonstructural protein 14 (NSP14), main protease (Mpro), and papain-like protease (PLpro) had the best docking scores with fluvoxamine. Consistent with the docking results, it was confirmed by molecular dynamics simulations that the NSP14 N7-MTase ((N7-guanine)-methyltransferase)-fluvoxamine, Mpro-fluvoxamine, and PLpro-fluvoxamine complexes are stable, with the lowest binding free energies of -105.1, -82.7, and - 38.5 kJ/mol, respectively. A number of hotspot residues involved in the interaction were also identified. These include Glu166, Asp187, His41, and Cys145 in Mpro, Gly163 and Arg166 in PLpro, and Glu302, Gly333, and Phe426 in NSP14, which could aid in the development of better antivirals against SARS-CoV-2.

On the Recognition of Natural Substrate CTP and Endogenous Inhibitor ddhCTP of SARS-CoV-2 RNA-Dependent RNA Polymerase: A Molecular Dynamics Study.

The novel coronavirus SARS-CoV-2 is the causative agent of the COVID-19 outbreak that is affecting the entire planet. As the pandemic is still spreading worldwide, with multiple mutations of the virus, it is of interest and of help to employ computational methods for identifying potential inhibitors of the enzymes responsible for viral replication. Attractive antiviral nucleotide analogue RNA-dependent RNA polymerase (RdRp) chain terminator inhibitors are investigated with this purpose. This study, based on molecular dynamics (MD) simulations, addresses the important aspects of the incorporation of an endogenously synthesized nucleoside triphosphate, ddhCTP, in comparison with the natural nucleobase cytidine triphosphate (CTP) in RdRp. The ddhCTP species is the product of the viperin antiviral protein as part of the innate immune response. The absence of the ribose 3'-OH in ddhCTP could have important implications in its inhibitory mechanism of RdRp. We built an in silico model of the RNA strand embedded in RdRp using experimental methods, starting from the cryo-electron microscopy structure and exploiting the information obtained by spectrometry on the RNA sequence. We determined that the model was stable during the MD simulation time. The obtained results provide deeper insights into the incorporation of nucleoside triphosphates, whose molecular mechanism by the RdRp active site still remains elusive.

Molecular dynamics simulations highlight the altered binding landscape at the spike-ACE2 interface between the Delta and Omicron variants compared to the SARS-CoV-2 original strain.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) B.1.1.529 variant (Omicron), represents a significant deviation in genetic makeup and function compared to previous variants. Following the BA.1 sublineage, the BA.2 and BA.3 Omicron subvariants became dominant, and currently the BA.4 and BA.5, which are quite distinct variants, have emerged. Using molecular dynamics simulations, we investigated the binding characteristics of the Delta and Omicron (BA.1) variants in comparison to wild-type (WT) at the interface of the spike protein receptor binding domain (RBD) and human angiotensin converting enzyme-2 (ACE2) ectodomain. The primary aim was to compare our molecular modelling systems with previously published observations, to determine the robustness of our approach for rapid prediction of emerging future variants. Delta and Omicron were found to bind to ACE2 with similar affinities (-39.4 and -43.3 kcal/mol, respectively) and stronger than WT (-33.5 kcal/mol). In line with previously published observations, the energy contributions of the non-mutated residues at the interface were largely retained between WT and the variants, with F456, F486, and Y489 having the strongest energy contributions to ACE2 binding. Further, residues N440K, Q498R, and N501Y were predicted to be energetically favourable in Omicron. In contrast to Omicron, which had the E484A and K417N mutations, intermolecular bonds were detected for the residue pairs E484:K31 and K417:D30 in WT and Delta, in accordance with previously published findings. Overall, our simplified molecular modelling approach represents a step towards predictive model systems for rapidly analysing arising variants of concern.

Characterization of Entry Pathways, Species-Specific Angiotensin-Converting Enzyme 2 Residues Determining Entry, and Antibody Neutralization Evasion of Omicron BA.1, BA.1.1, BA.2, and BA.3 Variants.

The SARS-CoV-2 Omicron variants were first detected in November 2021, and several Omicron lineages (BA.1, BA.2, BA.3, BA.4, and BA.5) have since rapidly emerged. Studies characterizing the mechanisms of Omicron variant infection and sensitivity to neutralizing antibodies induced upon vaccination are ongoing by several groups. In the present study, we used pseudoviruses to show that the transmembrane serine protease 2 (TMPRSS2) enhances infection of BA.1, BA.1.1, BA.2, and BA.3 Omicron variants to a lesser extent than ancestral D614G. We further show that Omicron variants have higher sensitivity to inhibition by soluble angiotensin-converting enzyme 2 (ACE2) and the endosomal inhibitor chloroquine compared to D614G. The Omicron variants also more efficiently used ACE2 receptors from 9 out of 10 animal species tested, and unlike the D614G variant, used mouse ACE2 due to the Q493R and Q498R spike substitutions. Finally, neutralization of the Omicron variants by antibodies induced by three doses of Pfizer/BNT162b2 mRNA vaccine was 7- to 8-fold less potent than the D614G. These results provide insights into the transmissibility and immune evasion capacity of the emerging Omicron variants to curb their ongoing spread. IMPORTANCE The ongoing emergence of SARS-CoV-2 Omicron variants with an extensive number of spike mutations poses a significant public health and zoonotic concern due to enhanced transmission fitness and escape from neutralizing antibodies. We studied three Omicron lineage variants (BA.1, BA.2, and BA.3) and found that transmembrane serine protease 2 has less influence on Omicron entry into cells than on D614G, and Omicron exhibits greater sensitivity to endosomal entry inhibition compared to D614G. In addition, Omicron displays more efficient usage of diverse animal species ACE2 receptors than D614G. Furthermore, due to Q493R/Q498R substitutions in spike, Omicron, but not D614G, can use the mouse ACE2 receptor. Finally, three doses of Pfizer/BNT162b2 mRNA vaccination elicit high neutralization titers against Omicron variants, although the neutralization titers are still 7- to 8-fold lower those that against D614G. These results may give insights into the transmissibility and immune evasion capacity of the emerging Omicron variants to curb their ongoing spread.

Intermediates in SARS-CoV-2 spike-mediated cell entry.

SARS-CoV-2 cell entry is completed after viral spike (S) protein-mediated membrane fusion between viral and host cell membranes. Stable prefusion and postfusion S structures have been resolved by cryo-electron microscopy and cryo-electron tomography, but the refolding intermediates on the fusion pathway are transient and have not been examined. We used an antiviral lipopeptide entry inhibitor to arrest S protein refolding and thereby capture intermediates as S proteins interact with hACE2 and fusion-activating proteases on cell-derived target membranes. Cryo-electron tomography imaged both extended and partially folded intermediate states of S2, as well as a novel late-stage conformation on the pathway to membrane fusion. The intermediates now identified in this dynamic S protein-directed fusion provide mechanistic insights that may guide the design of CoV entry inhibitors.

Shared N417-Dependent Epitope on the SARS-CoV-2 Omicron, Beta, and Delta Plus Variants.

As severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to evolve, several variants of concern (VOCs) have arisen which are defined by multiple mutations in their spike proteins. These VOCs have shown variable escape from antibody responses and have been shown to trigger qualitatively different antibody responses during infection. By studying plasma from individuals infected with either the original D614G, Beta, or Delta variants, we showed that the Beta and Delta variants elicit antibody responses that are overall more cross-reactive than those triggered by D614G. Patterns of cross-reactivity varied, and the Beta and Delta variants did not elicit cross-reactive responses to each other. However, Beta-elicited plasma was highly cross-reactive against Delta Plus (Delta+), which differs from Delta by a single K417N mutation in the receptor binding domain, suggesting that the plasma response targets the N417 residue. To probe this further, we isolated monoclonal antibodies from a Beta-infected individual with plasma responses against Beta, Delta+, and Omicron, which all possess the N417 residue. We isolated an N417-dependent antibody, 084-7D, which showed similar neutralization breadth to the plasma. The 084-7D MAb utilized the IGHV3-23*01 germ line gene and had somatic hypermutations similar to those of previously described public antibodies which target the 417 residue. Thus, we have identified a novel antibody which targets a shared epitope found on three distinct VOCs, enabling their cross-neutralization. Understanding antibodies targeting escape mutations, such as K417N, which repeatedly emerge through convergent evolution in SARS-CoV-2 variants, may aid in the development of next-generation antibody therapeutics and vaccines. IMPORTANCE The evolution of SARS-CoV-2 has resulted in variants of concern (VOCs) with distinct spike mutations conferring various immune escape profiles. These variable mutations also influence the cross-reactivity of the antibody response mounted by individuals infected with each of these variants. This study sought to understand the antibody responses elicited by different SARS-CoV-2 variants and to define shared epitopes. We show that Beta and Delta infections resulted in antibody responses that were more cross-reactive than the original D614G variant, but they had differing patterns of cross-reactivity. We further isolated an antibody from Beta infection which targeted the N417 site, enabling cross-neutralization of Beta, Delta+, and Omicron, all of which possess this residue. The discovery of antibodies which target escape mutations common to multiple variants highlights conserved epitopes to target in future vaccines and therapeutics.

The Glycan-Binding Trait of the Sarbecovirus Spike N-Terminal Domain Reveals an Evolutionary Footprint.

The spike protein on sarbecovirus virions contains two external, protruding domains: an N-terminal domain (NTD) with unclear function and a C-terminal domain (CTD) that binds the host receptor, allowing for viral entry and infection. While the CTD is well studied for therapeutic interventions, the role of the NTD is far less well understood for many coronaviruses. Here, we demonstrate that the spike NTD from SARS-CoV-2 and other sarbecoviruses binds to unidentified glycans in vitro similarly to other members of the Coronaviridae family. We also show that these spike NTD (S-NTD) proteins adhere to Calu3 cells, a human lung cell line, although the biological relevance of this is unclear. In contrast to what has been shown for Middle East respiratory syndrome coronavirus (MERS-CoV), which attaches sialic acids during cell entry, sialic acids present on Calu3 cells inhibited sarbecovirus infection. Therefore, while sarbecoviruses can interact with cell surface glycans similarly to other coronaviruses, their reliance on glycans for entry is different from that of other respiratory coronaviruses, suggesting sarbecoviruses and MERS-CoV have adapted to different cell types, tissues, or hosts during their divergent evolution. Our findings provide important clues for further exploring the biological functions of sarbecovirus glycan binding and adds to our growing understanding of the complex forces that shape coronavirus spike evolution. IMPORTANCE Spike N-terminal domains (S-NTD) of sarbecoviruses are highly diverse; however, their function remains largely understudied compared with the receptor-binding domains (RBD). Here, we show that sarbecovirus S-NTD can be phylogenetically clustered into five clades and exhibit various levels of glycan binding in vitro. We also show that, unlike some coronaviruses, including MERS-CoV, sialic acids present on the surface of Calu3, a human lung cell culture, inhibit SARS-CoV-2 and other sarbecoviruses. These results suggest that while glycan binding might be an ancestral trait conserved across different coronavirus families, the functional outcome during infection can vary, reflecting divergent viral evolution. Our results expand our knowledge on the biological functions of the S-NTD across diverse sarbecoviruses and provide insight on the evolutionary history of coronavirus spike.

Omicron Binding Mode: Contact Analysis and Dynamics of the Omicron Receptor-Binding Domain in Complex with ACE2.

On 26 November 2021, the WHO classified the Omicron variant of the SARS-CoV-2 virus (B.1.1.529 lineage) as a variant of concern (VOC) (COVID-19 Variant Data, Department of Health, 2022). The Omicron variant contains as many as 26 unique mutations of effects not yet determined (Venkatakrishnan, A., Open Science Framework, 2021). Out of its total of 34 Spike protein mutations, 15 are located on the receptor-binding domain (S-RBD) (Stanford Coronavirus Antiviral & Resistance Database, 2022) that directly contacts the angiotensin-converting enzyme 2 (ACE2) host receptor and is also a primary target for antibodies. Here, we studied the binding mode of the S-RBD domain of the Spike protein carrying the Omicron mutations and the globular domain of human ACE2 using molecular dynamics (MD) simulations. We identified new and key Omicron-specific interactions such as R493 (of mutation Q493R), which forms salt bridges both with E35 and D38 of ACE2, Y501 (N501Y), which forms an edge-to-face aromatic interaction with Y41, and Y505 (Y505H), which makes an H-bond with E37 and K353. The glycan chains of ACE2 also bind differently in the WT and Omicron variants in response to different charge distributions on the surface of Spike proteins. However, while the Omicron mutations considerably improve the overall electrostatic fit of the two interfaces, the total number of specific and favorable interactions between the two does not increase. The dynamics of the complexes are highly affected too, making the Omicron S-RBD:ACE2 complex more rigid; the two main interaction sites, Patches I and II, isolated in the WT complex, become connected in the Omicron complex through the alternating interaction of R493 and R498 with E35 and D38.

Reconstruction of the unbinding pathways of noncovalent SARS-CoV and SARS-CoV-2 3CLpro inhibitors using unbiased molecular dynamics simulations.

The main protease (3CLpro) is one of the essential components of the SARS-CoVs viral life cycle, which makes it an interesting target for overpowering these viruses. Although many covalent and noncovalent inhibitors have been designed to inhibit this molecular target, none have gained FDA approval as a drug. Because of the high rate of COVID-19 pandemic development, in addition to laboratory research, we require in silico methods to accelerate rational drug design. The unbinding pathways of two SARS-CoV and SARS-CoV-2 3CLpro noncovalent inhibitors with the PDB IDs: 3V3M, 4MDS, 6W63, 5RF7 were explored from a comparative perspective using unbiased molecular dynamics (UMD) simulations. We uncovered common weak points for selected inhibitors that could not interact significantly with a binding pocket at specific residues by all their fragments. So water molecules entered the free binding S regions and weakened protein-inhibitor fundamental interactions gradually. N142, G143, and H163 are the essential residues, which cause key protein-ligand interactions in the binding pocket. We believe that these results will help design new potent inhibitors against SARS-CoV-2.

Kinetic Characterization of SARS-CoV-2 nsp13 ATPase Activity and Discovery of Small-Molecule Inhibitors.

SARS-CoV-2 non-structural protein 13 (nsp13) is a highly conserved helicase and RNA 5'-triphosphatase. It uses the energy derived from the hydrolysis of nucleoside triphosphates for directional movement along the nucleic acids and promotes the unwinding of double-stranded nucleic acids. Nsp13 is essential for replication and propagation of all human and non-human coronaviruses. Combined with its defined nucleotide binding site and druggability, nsp13 is one of the most promising candidates for the development of pan-coronavirus therapeutics. Here, we report the development and optimization of bioluminescence assays for kinetic characterization of nsp13 ATPase activity in the presence and absence of single-stranded DNA. Screening of a library of 5000 small molecules in the presence of single-stranded DNA resulted in the discovery of six nsp13 small-molecule inhibitors with IC50 values ranging from 6 ± 0.5 to 50 ± 6 µM. In addition to providing validated methods for high-throughput screening of nsp13 in drug discovery campaigns, the reproducible screening hits we present here could potentially be chemistry starting points toward the development of more potent and selective nsp13 inhibitors, enabling the discovery of antiviral therapeutics.

Differential Interactions between Human ACE2 and Spike RBD of SARS-CoV-2 Variants of Concern.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the current coronavirus disease 2019 (COVID-19) pandemic. It is known that the receptor-binding domain (RBD) of the spike protein of SARS-CoV-2 interacts with the human angiotensin-converting enzyme 2 (ACE2) receptor, initiating the entry of SARS-CoV-2. Since its emergence, a number of SARS-CoV-2 variants have been reported, and the variants that show high infectivity are classified as variants of concern according to the United States Centers for Disease Control and Prevention. In this study, we performed both all-atom steered molecular dynamics (SMD) simulations and microscale thermophoresis (MST) experiments to characterize the binding interactions between ACE2 and RBD of all current variants of concern (Alpha, Beta, Gamma, and Delta) and two variants of interest (Epsilon and Kappa). We report that RBD of the Alpha (N501Y) variant requires the highest amount of force initially to be detached from ACE2 due to the N501Y mutation in addition to the role of N90-glycan, followed by Beta/Gamma (K417N/T, E484 K, and N501Y) or Delta (L452R and T478 K) variants. Among all variants investigated in this work, RBD of the Epsilon (L452R) variant is relatively easily detached from ACE2. Our results from both SMD simulations and MST experiments indicate what makes each variant more contagious in terms of RBD and ACE2 interactions. This study could shed light on developing new drugs to inhibit SARS-CoV-2 entry effectively.

Pathogen-sugar interactions revealed by universal saturation transfer analysis.

Many pathogens exploit host cell-surface glycans. However, precise analyses of glycan ligands binding with heavily modified pathogen proteins can be confounded by overlapping sugar signals and/or compounded with known experimental constraints. Universal saturation transfer analysis (uSTA) builds on existing nuclear magnetic resonance spectroscopy to provide an automated workflow for quantitating protein-ligand interactions. uSTA reveals that early-pandemic, B-origin-lineage severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike trimer binds sialoside sugars in an "end-on" manner. uSTA-guided modeling and a high-resolution cryo-electron microscopy structure implicate the spike N-terminal domain (NTD) and confirm end-on binding. This finding rationalizes the effect of NTD mutations that abolish sugar binding in SARS-CoV-2 variants of concern. Together with genetic variance analyses in early pandemic patient cohorts, this binding implicates a sialylated polylactosamine motif found on tetraantennary N-linked glycoproteins deep in the human lung as potentially relevant to virulence and/or zoonosis.

A So-Far Overlooked Secondary Conformation State in the Binding Mode of SARS-CoV-2 Spike Protein to Human ACE2 and Its Conversion Rate Are Crucial for Estimating Infectivity Efficacy of the Underlying Virus Variant.

Since its outbreak in 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread with high transmission efficiency across the world, putting health care as well as economic systems under pressure. During the course of the pandemic, the originally identified SARS-CoV-2 variant has been multiple times replaced by various mutant versions, which showed enhanced fitness due to increased infection and transmission rates. In order to find an explanation for why SARS-CoV-2 and its emerging mutated versions showed enhanced transmission efficiency compared with SARS-CoV (2002), an enhanced binding affinity of the spike protein to human angiotensin converting enzyme 2 (hACE2) has been proposed by crystal structure analysis and was identified in cell culture models. Kinetic analysis of the interaction of various spike protein constructs with hACE2 was considered to be best described by a Langmuir-based 1:1 stoichiometric interaction. However, we demonstrate in this report that the SARS-CoV-2 spike protein interaction with hACE2 is best described by a two-step interaction, which is defined by an initial binding event followed by a slower secondary rate transition that enhances the stability of the complex by a factor of ~190 (primary versus secondary state) with an overall equilibrium dissociation constant (KD) of 0.20 nM. In addition, we show that the secondary rate transition is not only present in SARS-CoV-2 wild type ("wt"; Wuhan strain) but also found in the B.1.1.7 variant, where its transition rate is 5-fold increased. IMPORTANCE The current SARS-CoV-2 pandemic is characterized by the high infectivity of SARS-CoV-2 and its derived variants of concern (VOCs). It has been widely assumed that the reason for its increased cell entry compared with SARS-CoV (2002) is due to alterations in the viral spike protein, where single amino acid residue substitutions can increase affinity for hACE2. So far, the interaction of a single unit of the CoV-2 spike protein has been described using the 1:1 Langmuir interaction kinetic. However, we demonstrate here that there is a secondary state binding step that may be essential for novel VOCs in order to further increase their infectivity. These findings are important for quantitatively understanding the infection process of SARS-CoV-2 and characterization of emerging SARS-CoV-2 variants of spike proteins. Thus, they provide a tool for predicting the potential infectivity of the respective viral variants based on secondary rate transition and secondary complex stability.